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Eigenstates in finite systems such as nuclei, atoms, atomic clusters and quantum dots with few 
excited particles are chaotic superpositions of shell model basis states. We study criterion for the 
equilibrium distribution of basis components (ergodicity, or Quantum Chaos), effects of level density 
variation and transition from the Breit-Wigner to the Gaussian shape of eigenstates and strength 
functions. In the model of n interacting particles distributed over m orbitals, the shape is given by 
the Breit-Wigner function with the width in the form of gaussian dependence on energy. 

PACS numbers: 05.45. +b, 3I.25.-v, 31.50.+W, 32.30.-r 



Introduction. Recently, based on chaotic structure 
of eigenstates, statistical approach has been developed in 
which allows to find distribution of occupation num- 
bers for single-particles states, expectation values of dif- 
ferent operators, transition amplitudes between chaotic 
states and degree of an enhancement of weak interactions. 
The main quantity in this approach is the form of the 
distribution of shell model basis components in chaotic 
eigenstates. In recent numerical studies of the Ce atom 
H , the s—d nuclear shell model ]Bj and random two-body 
interaction model [pl[| it was found that typical shape of 
exact eigenstates practically does not depend on a partic- 
ular many-body system and has a universal form which 
essentially depends on few parameters. 

In this paper, we develop a method for the description 
of the form of chaotic eigenstates and strength functions 
without diagonalization of huge many-body Hamiltonian 
matrices. 

We follow the standard way according to which the to- 
tal Hamiltonian for n Fermi-particles distributed over m 
single-particles states is written as a sum of two terms, 



H = H Q + V = e s ata s + 

s=l 



-^2V pqrs a+a+a r a s (I) 



Here the "unperturbed" Hamiltonian Hq incorporates an 
effect of the mean field, e s are the energies of single- 
particle states calculated in this field, af, a s are creation- 
annihilation operators, and V stands for the residual two- 
body interaction (the difference between the exact and 
mean field Hamiltonians). 

Exact ("compound") eigenstates |i) of the Hamiltonian 
H can be expressed in terms of simple shell-model basis 
states \k) (eigenstates of Hq) : 



\i) = 



(i) 



\k) ; \k) 



10) 



(2) 



These compound eigenstates \i) are formed by the resid- 
ual interaction V ; in complex systems they typically 



contain large number N pc 3> I of the so-called principal 

(i) 

components C k which fluctuate "randomly" as a func- 
tion of indices i and k. The shape of exact eigenstates 
is given by the "spreading function" F (in what follows, 
the F — function ), 



F, 



(0 - 



C, 



(i) 



(3) 



where E^ is the unperturbed energy and E^ is the per- 
turbed one. 

Equations for strength function and spreading 
width. For a weak interaction between particles the 
shape of chaotic eigenstates is known to be well described 
by the Breit-Wigner form. However, with an increase of 
the interaction strength, the average shape of the eigen- 
states (F— function) changes from the Breit-Wigner one 
to that close to the Gaussian with the exponential tails 
. In order to reduce the distortion effect due to non- 
constant density of states, in what follows we consider the 
so-called "strength function" (E) which is also known 
in literature as "local spectral density of states" (LDOS), 



P k (E)=F(E k ,E)p{E) 



(4) 



Note that the F— function gives the shapes of both exact 
eigenstates and the LDOS, depending on what is fixed, 
the total energy E^ or the unperturbed one, E k . 

The equations for P k (E) can be obtained on the base of 
the method presented in Let us choose some basis 

component |fc) and diagonalize the Hamiltonian matrix 
without this component. Thus, the problem is reduced 
to the interaction of this component with the exact eigen- 
states \i) described by the matrix elements V k i- 

The solution can be written by making use of the av- 
erage over a small energy interval A (see details in H), 



Pk(E) = 



1 



r k (E) 



2ir (E k +5 k -EY + {T k {E)/2f 



(5) 



1 



T k (E) ~2ii\V kl \ 2 p(E) 



E 



\Vki\ 2 (E-EV) 
(E - EWy + (A/2) 2 



(6) 



(7) 



It is easy to recognize in the energy shift S k the mod- 
ified second order correction to the unperturbed energy 
level. For the calculation of the shape of the eigenvector 
\i > one should substitute the exact energy E = E^' — 
Ei + 5i . The difference Si — S k is negligible if the interac- 
tion is not very strong. 

One should stress that the summation in the above 
equations are performed over exact states. We would like 
to express the result in terms of the basis states only, this 
allows us to solve the problem without diagonalization of 
the Hamiltonian matrix H kp . To do this, we express ex- 
act eigenstates \i > in terms of the basis components, 



\Vki 



\Hkp\ 2 +Y J C q >C p )H k P H qk 



(8) 



with Hkp standing for non-diagonal Hamiltonian matrix 
elements. If coefficients Cp/ q can be treated as random 
variables, the second term vanishes after averaging. Sub- 
stitution of Eq. (||) into Eqs.(^, fj]) gives 

T k (E) = 2nY,\Hkp\ 2 P P (E) = 



p^k 



2_j \ H kp\ 
P¥=k 



?p{E) 



{Ep + 8p-EY + {Tp(E)/2f 



(9) 



Sk — ^ \H kp \ 
p^k 



dE ( ^ 



P P (E^) 
E - Ed) 



\Hk P \ 2 {E — E p — S p ) 



P ^k 



'_{E-E p -S p ) 2 + {T p {E)/2f 



(10) 



where the integral is taken as the principal value. Last 
equality is valid in the approximation of slow variation 
of T P (E) and 5 p . The equations for T k (E) and S k allow 
to calculate the strength function (||) from the unper- 
turbed energy spectrum and matrix elements of the total 
Hamiltonian H. 

Criterion of statistical equilibrium. Let us discuss 
the conditions under which self-consistent solution of the 
equations (j^, |^, ^) exists. There are four important pa- 
rameters in this problem: the spreading width of a basis 
component T, the effective band width a of the Hamilto- 
nian matrix H pq , the interval between the many - body 
energy levels D — p~ x and the interval df — pj 1 between 
the final basis states \p > which can be connected with a 
particular basis component \k > by the two-body inter- 
action. The ratio D/df is exponentially small M since all 



the basis states \p > which differ from \k > by position of 
more than two particles, have zero matrix elements H kp 
and do not contribute to the Eqs. (||, |i~0| ). 

First, let us consider Eqs.(||, [l(]) for a strong enough 
interaction, T » df. In this case the number of effec- 
tively large terms in the sums is large, Nf ~ T/d/, and 
fluctuations of T are small, 5T ~ Y/ yjNf. Therefore, Eq. 
(0) can be written as 



T k (E) ~ 2ir\H kp \ 2 pf(E) 



(11) 



< S p > can 



where E = E — 5. The energy shift S 
be neglected in the case of T << a. In order to per- 
form the summation over p, we assumed that r(Z?) and 
Pf(E) vary slowly within the energy interval of the size 
r. Thus, in order to have large number of final states 
Nf ~ 2irH 2 p /d 2 f and statistical equilibrium (small fluc- 
tuations of L), one needs H kp » df. In this case chaotic 
components of exact eigenfunctions in the unperturbed 
many-particle basis ergodically fill the whole energy shell 
of the width T, with Gaussian fluctuations of the coef- 
ficients cjp with the variance given by the F— function 
© (see also §1). 

With the decrease of the ratio H kp /df the fluctua- 
tions of r increase and for H kp < df the smooth self- 
consistent solution of Eqs.([)]) disappears. Indeed, in this 
case T p in the denominator can be neglected and the sum 
in (Q) is dominated by one term with a minimal energy 
E — E p ~ df. Therefore, for a typical basis state \k > 
formally one gets T k ~ Tp(H kp /df) 2 « T p . This con- 
tradicts to the equilibrium condition according to which 
all components are "equal" ( T k ~ T p ). 

One should stress that the absence of a smooth solution 
for the shape of the eigenstates and the strength function 
does not mean that the number of principal components 
in exact eigenstates is small. However, the distribution 
of the components is not ergodic: there are many "holes" 
inside exact eigenstates which occupy the energy shell of 

the width 2n\H kp \ 2 p f (E) (see @,@). In such a situa- 
tion, very large (non-Gaussian) fluctuations of are 
typical. 

It is important that ensemble averaging in this prob- 
lem is not equivalent to the energy average for a specific 
Hamiltonian matrix. For example, the average over the 
single-particle spectrum leads to variation of energy de- 
nominators in (j^) and can fill the holes in the F— function 
even for L < df. 

Transition from the Breit-Wigner type to the 
Gaussian-like strength function. In principal, the 
set of equations (|], [n]) for the shape of the strength 
function P k (E) defined by Eq.(|), can be solved numer- 
ically having the unperturbed many-body spectrum and 
matrix elements H kp of the total Hamiltonian. However, 
for relatively large number of particles (practically, for 
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n > 4), one can find an approximate analytical solution 
of the problem. 

First, we note that the spreading width T(E) in the 
expression (|^) for the strength function can be a strong 
function of excitation energy E due to the variation of 
the density of the final states pf(E)) — (df)~ x in Eq. 
(pi]). It is well-known that at small E the basis compo- 
nent with one excited particle has T(E) cx E 2 /do where 
<io is the interval between single-particle energy levels. 
For typical case of n* ~ (E/do) 1 ^ 2 excited particles the 
spreading width can be estimated as T(E) oc ~ 
(d )^ 1 (E/d ) 3 / 2 |§- Below we show that at higher en- 
ergies far from the ground state, the energy dependence 
of Pf{E) and T(E) can be quite close to the Gaussian. 
Note that the Gaussian form typically occurs in "statis- 
tical spectroscopy" when neglecting the mean field 
term in Eq. (]]]). 

In the model (|l]) the density df is defined by transi- 
tions between those basis states which differ by the po- 
sition of one or two particles only, therefore, Pf(E) — 

pf(E)+pf(E). 

(2) 

Let us estimate the density p^ determined by the en- 

(2) 

ergy difference ui pj ! between the states \p > and \k > 
which differ by the position of two particles, 




+ {Vau+Vp v -V rv -V 5u ) + V aP -V r s (12) 

Here the summation is taken over n — 2 occupied orbitals 
and V av are the diagonal matrix elements of the residual 
interaction between the particles located at the orbitals 
a and v. The matrix elements of residual interaction are 
assumed to be random with the zero mean. 

For large number of fluctuating terms in the Eq.jl^) 
the distribution of to is close to the gaussian form. 
Strictly speaking, this is correct if the contribution of 
An — 6 interaction terms to the frequency Eq.(|l2|) is 
strong. However, even four single-particle energy terms 
give the distribution which is close to the Gaussian. The 
same conclusion is reasonable also for single-particle tran- 
sition density p^\ thus, the general expression reads as 

p^\E) g K(2na>)- 1/2 e X p(- ^ = W? )■ (13) 

The normalization parameter K stands for the number of 
one or two-particle transitions, K — K\ = n(m — n) and 
K = K2 = n(n — l)(m — n)(m — n — l)/4 correspondingly 

From Eq.(|12[) the estimate for the average frequency 
of two-particle transitions reads as a/ 2 ) = 2(e~p~ — e^) w 
2m/(m — n)(e — Ek/n) where Ck = Ek/n is the average 
single-particle energy in the basis state fc > containing 
n particles, e is the single-particle energy averaged over 



all m orbitals. Average energy of the empty orbitals e p 
can be found from the relation me = e^n + e~p~(m — n) . 

The variance of pf'(E) for two-particle transitions is 
equal to erf = 2a 2 + 2a 2 + (An - 6)V 2 « 2a 2 + (An - 6)V 2 
where a 2 is the variance of single-particle spectrum, and 
V is the variance of non-diagonal matrix elements of the 
two-body residual interaction. Note that in the case of 
n « m for low- lying states the variance of the occupied 
orbital energies a\ is small and the variance of empty 
orbital energies is a 2 ~ a 2 . 

Similar, the density pQ' is also approximated by the 

Eq.(|3|), with K = K t , w m/(m-n)(e- Ek/n) and 
<A = °l + 4 + 2(n - l)V 2 w a 2 + 2(n - l)V 2 . 

Thus, the width V(E) is given by the following expres- 
sion, T = 27r[(n - l)V 2 p ( p + V 2 pf ] }. Since for single- 
particle transitions the summation in H^ p = ^ v V al ^^~ (L , 
is performed over occupied orbitals, the factor n — 1 
appears in the above relation. Typically, the ratio 
Kz/Kn — = (m — n — l)/4 is larger than 1, there- 

fore, the two-particle transitions dominate. In this case 
we can neglect the differences in ZZJ and a for two-particle 
and one-particle transitions and combine two terms into 
one. As a result, the spreading width is described by the 
simple Gaussian formula 

ri(E) , 2 . (aE) ^ rap {_(i^p£} 

(14) 

where E — E — 5. Here (AE) 2 is the variance of the 
strength function which can be defined through its aver- 
age value §j§, 

(AE) 2 k = H l P = V 2 n{n- l)(m- n)(3 + m- n)/A 

and u>k and Cfe are close to that for the two-particle tran- 
sitions. The maximum of Pf(E) and r(_E) is shifted by 
\ulk\ towards the center of the spectrum compared to the 
maximum of Breit-Wigner function. This leads to some 
distortion of the strength function Eq.(||) and the shape 
of the eigenstates, which is especially large at the bottom 
of the spectrum. 

Thus, we have demonstrated that if the interaction is 
small (r <C Tfe), the strength function has the Breit- 
Wigner shape with the broad gaussian envelope originat- 
ing from Tk(E) in the numerator of Eq. (|5|). It is easy 
to check that this envelope is, indeed, needed in order to 
provide the correct value (AE) 2 for the second moment of 
the strength function (note, that the Breit-Wigner shape 
has infinite second moment which is unphysical). 

When the interaction V increases one needs to take 
into account one more contribution to the broadening 
of the shape of F(E). It is given by the width of the 
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strength function P p {E) in Eq.(||) (it was neglected in 
Eq.(|i"l|)). Taking into account this width we can give an 
estimate a\ ~ a%+r 2 . With further increase of interac- 
tion, where the shape of P p (E) is close to the Gaussian, 
we have a\ ~ o\ + (AE) 2 k . 

Direct numerical study of the model (|l|) with n = 6 
Fermi-particles and m — 13 orbitals shows that the above 
analytical expressions give quite a good description of 
the shape of the strength function Pk(E) as well as of 
the energy dependence of the spreading width. The un- 
perturbed single-particle spectrum has been chosen at 
random, with do = I and e s ~ dgs. The size of the 
Hamiltonian matrix is N — = 1716 and we specify 
the unperturbed state \io > with io = 280. 



(numerically) T(E),P k {E) and F(E,E k ), using T from 
Eq.(Q with a 2 ~ a 2 + {&-E) 2 k as the zero approxima- 
tion in the right-hand side of Eqs.(||, |To| ). 
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FIG. 1. The F— function (|3j) in the basis representation. 
Broken line is the result of numerical diagonalization of the 
Hamiltonian matrix; to reduce the fluctuations, we take the 
average over 10 matrices Hit with different two-body random 
interaction with V 2 ~ 0.15. Two smooth curves correspond 
to the computation of Eq.(^|) with F k (E) given by Eqs.fl|l(]) 
and by Eq. (|l4|) ; they practically coincide. The inset shows 
the dependence Fk(E) itself; full curve is the expression ([li]), 
the dashed curve is the computation from Eqs.(g, |Io| ). 

In the estimates above we assumed that T << a where 
a is the effective energy band width of the Hamiltonian 
matrix, see Eq.(jl4]). When r ~ er, the (Gaussian) varia- 
tion of T(E) in the numerator of the strength function in 
Eq.(||) becomes as important as the variation of the Breit- 
Wigner energy denominator (E - E k ) 2 + (T/2) 2 . At this 
point, r « cr, the transition from the Breit-Wigner type 
to Gaussian-type shape of the eigenstates takes place. 
We still can use Eqs.(||, |l^), and in order to calculate 
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